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Abstract 

We use the c = 1 matrix model to compute the potential energy V{C) for (the zero 
mode of) the RR scalar in two-dimensional type OB string theory. The potential is 
induced by turning on a background RR flux, which in the matrix model corresponds 
to unequal Fermi levels for the two types of fermions. Perturbatively, this leads to 
a linear runaway potential, but non-perturbative effects stabilize the potential, and 
we find the exact expression V{C) = ^ J da arccos[cos(C)/vT+~e~^™]. We also 
compute the finite-temperature partition function of the OB theory in the presence of 
flux. The perturbative expansion is T-dual to the analogous result in type OA theory, 
but non-perturbative effects (which depend on C) do not respect naive R 1/R 
duality. The model can also be used to study scattering amplitudes in background RR 
fluxes. 
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1 Introduction 



The stabilization of moduli is an important problem on the way towards constructing 
realistic models in string/M-theory. It has been appreciated in recent years that a 
simple way of stabilizing moduli is to turn on fluxes in the compactification manifold. 
For example, in constructions based on the type JIB string or F-theory, turning on 
fluxes induces a superpotential which generically fixes all complex structure moduli pp. 
However, the Kahler moduli and their associated RR partners remain untouched in this 
procedure, and stabilizing them seems to require non-perturbative effects More 
precisely, all potentials that can be induced at the perturbative level (by breaking 
super symmetry) are runaway potentials, without any stable stationary point. But after 
including non-perturbative effects [Sj, the potential can have (meta-)stable minima, of 
great cosmological interest [2]. 

In this paper, we study the potential for a particularly simple "modulus", the 
zero mode of the Ramond-Ramond (RR) scalar C of type OB "non-critical" string 
theory in two dimensions. This constant part of C is related to a shift symmetry 
that is analogous to the gauge symmetries of p-form fields in higher-dimensional string 
theories. But while in higher dimensions, potentials for the periods of these p-form 
gauge fields are sometimes hard to compute due to the lack of a fully non-perturbative 
definition of string/M-theory, the advantage of the two-dimensional setup is that there 
exists an accessible non-perturbative definition of the theory, the so-called c = 1 matrix 
model [SlEl- We will therefore be able to compute everything exactly. 

As we will see, the mechanism for generating the potential in the two-dimensional 
OB context is essentially similar to the situation in higher dimensions. One can induce 
a potential for C by turning on a RR background flux. Perturbatively, this potential 
would exhibit runaway, but this behavior is stabilized by non-perturbative effects. The 
exact expression is (see figure El on page Uni for a graph of this function) 



where /i is the tachyon background (worldsheet cosmological constant) determining the 
string coupling, and Q is the appropriately normalized RR flux. The crucial step in 
understanding V{C) is the identification, in the finite- matrix model, of the RR flux 
and the zero mode of C. Briefly, while n corresponds to the distance of the Fermi level 
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for the fermionic eigenvalues from the top of the inverted harmonic oscillator potential 
of the double-scaled matrix model, Q corresponds to the difference of the Fermi levels 
for the even and odd modes of the inverted harmonic oscillator^. Moreover, C is 
identified with the Z2-breaking boundary condition in the asymptotic regions of the 
potential. We will give a careful derivation of these identifications in the upcoming 
section |21 

In section El we then turn to the computation of the potential energy V{C) as a 
function of the string coupling and background flux. Since the finite piece vanishes 
at Q = 0, we will for completeness also compute the sub leading contribution that 
vanishes when we take the double-scaling limit. 

We then give two other applications of our identifications of section |21 In section 
m we compute the finite-temperature partition function for non-zero background flux 
Q, and compare with the analogous results from the type OA matrix model. We find 
that once the OA flux is continued to imaginary values, the perturbative expansions of 
the two results are T-dual to each other. The non-perturbative parts of the partition 
functions, however, are not mapped to each other under naive — > l/i? duality. We 
speculate on the meaning of this result. On the other hand, our results are invariant 
under "S-duality" , which involves /i — yU and dualization oi C ^ C. In addition, we 
find a novel duality which exchanges NS and RR background fi Q. In the matrix 
model, these dualities are simple particle/hole dualities for the even and odd modes 
of the potential. Some of the relevant computations, such as the dependence of the 
partition function on C, and the specialization of the results to the "self-dual" radius 
R = 1, can be found in the appendix. In sectional we study the bosonization of the 
matrix model eigenvalues when the Fermi levels for even and odd modes are not the 
same, and compute some S-matrix elements, following [7||H]. We conclude in section El 
with a list of open problems. 

For reviews of matrix models useful for our present discussion, see [11110111111121 • 
More recent matrix model literature includes [1311111131111111111111112011^11221121123 

i2nii2zii2Hii2niEni- 

^This is the non-perturbative definition. Perturbatively, it can be reduced to the difference of the 
Fermi levels on the two sides of the potential, or for left and right moving eigenvalues, depending on 
the sign of fi. This identification was proposed in (6). 
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2 The OB matrix model and its deformations 

The matrix model which in the double-scahng hmit describes two-dimensional M = 1 
supergravity coupled to c = 1 matter with type OB GSO projection [SHEl is essentially 
similar to the conventional or "old-fashioned" matrix model describing two-dimensional 
bosonic gravity coupled to c = 1 matter. The definition begins with the quantum me- 
chanics of an X hermitian matrix M with potential V{M) which has a quadratic 
maximum at M = 0. In the singlet sector, this quantum mechanics can be reduced 
to the dynamics the matrix' eigenvalues A behaving as free fermions moving in the 
potential V{X). The continuum limit involves taking to infinity, sending Planck's 
constant to zero, and adjusting the various parameters such that the Fermi level ap- 
proaches the top of the potential to order h. The only parameters surviving the limit 
are the Fermi level /i and the curvature of the potential at the maximum. The bosonic 
and the supersymmetric model differ slightly but importantly in both. 

In the "old-fashioned" way of thinking [^T1l32p33j , the expansion of the path- integral 
for finite A^ captures the discretization of the two-dimensional string worldsheet. The 
double-scaling limit is taken in order to tune the model to criticality. In the "reloaded" 
form Q , one thinks of the quantum mechanics as arising holographically as the world- 
volume theory on a stack of a large number of unstable DO-branes. In this context, the 
curvature of the potential at the maximum is identified with the mass of the open string 
tachyon. It is m? = — l/a' in the bosonic case and = —l/2a' in the supersymmetric 
case. 

Also, in the old days of the bosonic model, fermions were only filled on one side 
of the maximum, which can only make sense in the perturbative h expansion. For 
the supersymmetric version, we are now instructed to fill both sides of the potential, 
and to interpret this non-perturbatively stable model as the two-dimensional linear 
dilaton background of the type OB string. As in the bosonic case, collective excitations 
of the fermionic eigenvalues are identified with perturbative string theory degrees of 
freedom. In the bosonic case, there is only the center of mass of the string, which in 
two dimensions is a massless tachyon T. In the supersymmetric OB model, one has in 
addition a massless RR scalar C. According to [HlEl; even fiuctuations of the Fermi 
sea correspond to T, and odd fiuctuations correspond to C. In other words, the Z2 
parity symmetry of the matrix model M — > —M is identified with (—1)'^^ in the string 
theory picture. 
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In this paper, we will be interested in studying aspects of the RR sector of this 
OB matrix model. As we have mentioned in the introduction, the RR scalar shares 
many interesting features with the RR p-form fields familiar from higher-dimensional 
string theories. For type theories in general, the spacetime effective action contains 
to lowest order terms of the form [311IH3EniEZlin] 

l.Jf{T)dCA*dC (2) 

that couple the RR forms to the closed string tachyon T. Constraints on the functions 
f(T) can be obtained from the requirement of T-duality jSZj, and suggest taking f(T) = 
e~^"^, as we will do, following jH]. In two dimensions, there is only one RR form, which is 
a (non-chiral) scalar boson. As usual, the zero mode of C does not appear in the action 
and we have a one-parameter shift symmetry. Similarly to higher dimensions, this 
perturbative gauge symmetry is violated non-perturbatively by instanton effects. Such 
instantons should generically generate a potential for C, and by analogy with the Yang- 
Mills 6'-parameter, one might expect that this potential has to lowest order the form 
g-i/9a cosC. We will see that the details are in fact somewhat different in the present 
case. 

The natural field strength associated with C is F = e~^dC, satisfying the Bianchi 
identity and equation of motion 

d(e'^F) =0, 

(3) 

d{e-^ *F) =0. 

Recall that to lowest order, the tachyon background is given by 

T(0,t)=T(0)=/ie^ (4) 

where (j) is the spacelike Liouville direction supporting the linear dilaton. In this back- 
ground the linearized equations of motion for F admit the two linearly independent 
solutions 0, 

F = e-^dt , (5) 
F=e^d(t), (6) 

to which we shall refer to as the electric and magnetic solution, respectively. Of course, 
since C is the middle- dimensional form in the two-dimensional OB theory, we could 
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equally well have chosen to describe this using the dual form C related to C by e^^dC = 
F = *F = *dC, so that the action for C is (see [TTJ for some details on electric- 
magnetic dualities in type theories) 



We also note that depending on the sign of /i, one of the two solutions (0) decays 
at infinity ^ oo, which as we recall is the asymptotic region where the matrix model 
lives. For clarity, we will primarily discuss the case > 0, by which we mean that the 
Fermi level is below the top of the potential. In this case, only the electric solution 
(0) is physically relevant.^ But as is apparent from our discussion, one in facts expects 
that the theory enjoys an exact S-like duality /i —fi, C ^ C . We will confirm that 
after a proper identification of C and C, all our results are invariant under this duality. 

In it was proposed to identify the generator of the shift symmetry of C with 
the perturbatively (for yU > 0) well-defined difference in the number of fermions on 
the right and on the left of the maximum of the potential in the matrix model. One 
can justify this identification from the fact that instantons in spacetime are related to 
the tunneling of eigenvalues through the potential barrier in the matrix model. The 
instantons carry C-field charge and the tunneling changes the number of fermions on 
the left and right by one unit. Thus, a constant shift C ^ C + 2a shifts the phase of 
the wavefunction of fermions on the left of the potential by —a and the phase of the 
wavefunction of fermions on the right by +a. 

Another identification provided in [U] is that the decaying solution Q is related to 
the difference of the Fermi level for fermions on the left and the right of the potential, 
which perturbatively is again a well-defined notion. But since we are now allowing 
eigenvalue tunneling, one is naturally led to wonder about the non-perturbative stabil- 
ity of this flux background. From the outset, we notice that the total flux associated 
with the solution (0), Ft d(f) , as well as the energy, / F^ d(j), diverge as the volume 
of space. One the other hand, an instanton carries a quantized (electric) charge, and a 
single tunneling can only change the flux by a finite amount. Moreover, in the space- 
time picture the instanton corresponds to a D-brane which is localized in the Liouville 
direction at infinite 0, and has no zero mode. It is therefore not clear whether instan- 

^Our conventions are slightly different from the ones used in j^j, where the magnetic solution is 
physical when the Fermi level is below the top of the potential. In other words, the field C of [H| is 
our magnetic variable C. 




(7) 
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Figure 1: Perturbative energy levels in the inverted harmonic oscillator potential (finite N). 
For symmetric boundary conditions (C = 0), the spectrum is doubly degenerate (left panel). 
Changing C shifts the phase of fermions on the left and right of the potential. This raises 
the levels on the left and lowers them on the right (the right panel is for C = 1, where the 
periodicity of C is 27r). The energy shifts cancel out in the sum and C is perturbatively an 
exact symmetry in the double-scaling limit. 

tons are actually able to destroy the flux completely. We will show in this paper that 
while the instantons do affect the identification of the RR sector in the matrix model, 
the flux background (0) is in fact perfectly stable. 

We begin in the perturbative picture and consider the energy levels of the matrix 
eigenvalues close to the top of the potential, see Fig. [U These energy levels are deter- 
mined semiclassically by a systematic WKB expansion. To lowest order, we have the 
Bohr-Sommerfeld quantization condition for the n-th energy level en, 

r-A, 



7rnh= j y2(e„ - V{\))d\, 



where we integrate up to the classical turning point satisfying ^(A*) = e„. The 
boundary condition at infinity depends on the details of the potential in this non- 
universal region, but can at most contribute an overall phase to the WKB approxima- 
tion. If the boundary condition is Z2 symmetric, the energy levels on the left and right 
will match, ei^n = ^R,n , as shown on the left in Fig. ^ We also record the familiar 
expression for the asymptotic density of states p, which is obtained from (jH)) upon 
approximating V^(A) quadratically around with V"{0) = —1, 



P 



dn 



dt 



In 2TThn 

(9) 



2Tih 

Shifting the phase of the wavefunction by a corresponds to adding ha on the RHS of 
(jHl). In the convenient normalization a ~ C/2, a shift of C by 27r maps the n-th energy 
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Figure 2: Perturbatively, we can fill the energy levels asymmetrically on the left and right 
of the potential, but keep the difference of Fermi levels finite in the double-scaling limit. This 
induces a linear potential for C. 



level on the right to the (n — l)st and the n-th energy level on the left to the {n + l)st. 
More precisely, we find from 

dn 1 , , 

that shifting C changes e^/ji by 

deL/R _ , 1 



'11^ 



dC 2np 

The two contributions cancel in the sum, so that if the energy levels are filled to the 
same Fermi level on the left and on the right, shifting C does not cost any energy, even 
before taking the double-scaling limit. In the double-scaling limit, which involves p 
diverging as | ln^|, a finite difference of fermion number or a finite shift of C does not 
change the renormalized Fermi level /z, which is measured from the top of the potential 
in units of h. Moreover, since we neglect tunneling, the scattering amplitudes do not 
feel any mismatch of energy levels on the left and right of the potential. We conclude 
that C is perturbatively an exact symmetry in the double-scaling limit. 

These considerations suggest a simple way of inducing a non-vanishing potential 
for C. The idea is to fill the eigenvalues asymmetrically, on the left up to /xl, and 
on the right up to fiR, and to keep the difference p^i — I^r finite in the double-scahng 
limit, see Fig. |21 We will parameterize jJii = 1^ + Q and = n — Q. In the double- 
scaling limit, non-zero Q makes two contributions to the energy. There is first of all 
an infinite contribution that arises because we have lifted an infinite number (in the 
limit) of fermions by a finite amount. This contribution will diverge as usual as | ln^|, 
which is identified with the volume of space in the spacetime picture. But there is also 
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a finite contribution, which depends on C. Indeed, upon integrating we find the 
contribution 

V{C) = jdep {en{C) - e«(0)) = . (12) 

(The derivation is similar for the other sign of Q.) We note that this term is non-zero 
but finite (not infinite) in the double-scahng hmit. 

One may wonder what kind of term in the spacetime action could capture a con- 
tribution of the form (fT^ . Following the suggestion made in |Hj, we assume that Q is 
related to the fiux of the RR one-form field strength, which in our conventions is the 
electric solution (0), F = e~^dC oc Qe~^dt. We can then write (fT^ as 

V{C) oc (CFi(+oo) - CFt(-oo)) = J d<f> d^{CFt) , (13) 

which corresponds to a total derivative term in the action of the form 

I d{CF) . (14) 

This term (fT^ is reminiscent of Chern-Simons terms familiar from higher-dimensional 
string theories. It would be interesting to see whether one can check the presence of 
such a term from the spacetime or Liouville point of view. The potential (P) we will 
compute can be viewed as the non-perturbative completion of the expression (|THl . 

As we have mentioned, the potential energy ()12j) is linear in C and does not have 
any stationary point. In addition, the fiux background in Fig. |21 appears to be unstable 
to eigenvalue tunneling. We will now give a qualitative description of how both C and 
the fiux are stabilized once we include non-perturbative effects, and will present the 
quantitative calculations in the next section. 

The picture is quite simple (see Fig. E)). If we allow for tunneling, the energy levels 
on the left and right of the potential mix, and the true eigenstates of the inverted 
harmonic oscillator are linear combinations of the two. For symmetric boundary con- 
ditions (C = 0), the eigenstates split into even and odd under A —A, with splitting 
(at finite A^) of order e"'^", where a is the distance from the top of the potential in units 
of h. Each eigenfunction is equally supported on the left and right of the potential. 
When C is non-zero, the perturbative energy levels on the left and right do not match 
(see Fig.lH) and mixing is suppressed. Therefore, the eigenf unctions on each side of the 
potential essentially keep their identity, with small admixture from the other side. The 
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Figure 3: Energy eigenvalues of tiie inverted harmonic oscillator potential (at finite N) as 
a function of the boundary condition asymmetry C. 

eigenvalues depend linearly on C, as in eq. (fTT|). When C reaches vr, the perturbative 
energy levels would cross, leading to eigenvalue repulsion. The eigenfunctions are again 
equally supported on both sides of the potential. From C = n to C = 2n, the picture is 
the same with inverted sign, and the full periodicity is C ^ C + 2tt. This dependence 
of the eigenvalues on C is depicted in Fig. |21 

Even though for C 7^ 0, the eigenstates of the inverted harmonic oscillator are 
strictly speaking neither even nor odd, we will simply refer to the states by the parity 
they have once adiabatically continued to C = 0. An equivalent notion is the parity 
of the wavefunction under the combined reflection A —A, C —C. But one has to 
keep in mind that except for small regions of order e"'^"' around C = and tt, these 
wavefunctions are actually well supported on one side of the potential. 

Putting everything together, we are now led to the identification of the electric 
RR flux background Q with the difference of the Fermi levels for the even and odd 
eigenmodes of the inverted harmonic oscillator. Perturbatively, this definition reduces 
to the earlier one given in [0]. To be specific, we will put the Fermi level for the even 
modes at fi+ = fi + Q and for the odd modes at /x_ = /i — Q, and we assume for the 
time being that both yU+ and yU_ are positive. The RR background thus corresponds 
to a natural deformation of the OB matrix model, and is perfectly stable. Moreover, 
once Q is non-zero, we induce a potential energy V{C) for the zero mode of C, which 
is obtained by summing up the dependence of the unpaired even/odd eigenvalues from 
/i_ to Before we turn to the computation of V{C), we will give an independent 
heuristic argument for this identification of the RR flux in the matrix model. 

According to IHIHl, we can think of the matrix model holographically as the world- 
volume theory on a stack of N unstable DO-branes of type OB string theory. In this 
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context, the matrix M can be identified with the open string tachyon T of the DO- 
branes, and the inverted harmonic oscillator is the maximum of the tachyon potential 
V(T). As is by now quite well established, open string tachyons on unstable D-branes 
also couple to RR fields via terms of the form pEj 

jTTdJ^{T)AC. (15) 

The purpose of this coupling is to give the correct RR charge to a tachyonic kink 
interpolating between the two minima of V, with the BPS condition AjF = / V relating 
charge and mass. The couplings V and are not known in general — in particular, they 
might depend on the closed string tachyon in type OB theory. However, all we need 
to retain presently is the behavior near the maximum of V. Then, since the constant 
part of does not matter, all we need to know is that while V is generically an even 
function under T —T, T is generically an odd function. 

Eq. (|T5|l is also relevant for describing the behavior of unstable D-branes in an 
external field. Upon integration by parts, we obtain |3] 

TrJ^(T)AF (16) 

where F = dC is the field strength associated with C, and we note again that is odd 
under T — > — T. 

All these facts strongly suggest that turning on RR flux should correspond holo- 
graphically to a deformation of the matrix model which is odd under M —M (or 
equivalently, A —A). In second quantized language, we are led to the Hamiltonian 

H = jd\ [^^^^ - y ^^(A)^(A) + + gvi;t(A)vi;(_A)] . (17) 

The equations of motion that follow from (fTTj) for the even and odd parts of \& imme- 
diately imply the identification of even and odd Fermi levels /i± = /i ± Q that we have 
claimed above. 

Finally, we would like to make a minor comment concerning the nature of the flux 
appearing in (fT3j) and (fT^ . In general, type string theories which are constructed 
with M = 2 worldsheet supersymmetry contain twice as many RR fields as their type 
II cousins. The middle-dimensional form has both electric and magnetic degrees of 
freedom. This doubling arises because the RR ground states of the NSR string with 
both even and odd worldsheet fermion number contribute. In contrast, the present two- 
dimensional type OB theory actually only enjoys M = 1 worldsheet supersymmetry. 
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and the corresponding type IIB model is dead (but see [23j). (One is tempted to call 
this the |B string.) It is no surprise, therefore, to find that there is only one kind of 
allowed RR background, which with our choice of variables is the electric background 
©• 

An analogous discussion holds for D-branes. Type OB theories with M = 2 world- 
sheet supersymmetry have two kinds of D-branes, conventionally called Dp+ and 
Dp—, each of which is charged (electrically) under the appropriate combination of 
RR fields jnHESllSE] • The middle-dimensional form has one electrically and one mag- 
netically charged brane. In our case, we have a reduction also of the number of branes 
because we only have AT = 1 worldsheet supersymmetry [H]. 

To complete the story, we therefore have to justify that the RR form that couples to 
the tachyon on our unstable DO-branes is indeed the electric variable, and not possibly 
the magnetic one, C. Here, we can use the results of ref. [HHIj in which Sen's brane 
descent relations where studied for type theories. According to jSHj there are two 
completely independent descent charts, one for the + and one for the — branes. If 
we assume that the electric fiux is the difference of Fermi levels, the instanton which 
corresponds to tunneling of eigenvalues is charged electrically under C. Therefore, 
since the instanton descends from the Euclidean kink on the unstable DO-brane, we 
conclude that the unstable DO-brane must indeed couple to the electric variable C, as 
we had needed for consistency. 

3 Computation of the potential 

In our notation, we will mostly follow the usual conventions of [HUll40| l^ . except that 
we will call Planck's constant h instead of reserving (3 for the inverse temperature 
appearing in section |3J In the double-scaling limit, we take /i — > 0, keeping various 
other quantities fixed. The Schrodinger equation for the matrix model eigenvalue A is 
given by 

{-^dl + V{\)-e)i;{eA) = ^ (18) 

where V^(A) is the matrix model potential. In the double-scaling limit, only the 
quadratic behavior near the maximum is relevant. We will make the convenient 
choice jlOj 

V{\) = \{l-\'), (19) 
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with infinite walls (Diriclilet boundary conditions on -0) near A = ±1. For comparison 
with we note that we have chosen units in which a' = 1/2. 
It is convenient to introduce rescaled variables, 

fh 1 

A = Y —X and e = - — ha , (20) 
in terms of which eq. (|T8|l becomes 

2 

[^l + ^-a)^P{a,x)=0. (21) 

This equation is solved by the parabolic cylinder functions called W{a,x), W{a, —x) 
in the conventions of |46j . 

The zero mode of C can easily be implemented by asymmetrically fine-tuning the 
positions of the walls in the limit h ^ 0. Explicitly, we impose Dirichlet boundary 
conditions ip{a, xr) = ip^a, xl) = 0, where 

(This corresponds to A/, = — 1 + hC /2 and A/j = 1 + hC/2.) One can easily check 
that shifting the walls in this way shifts the perturbative energy levels as in (jllj) . The 
eigenvalue equation we have to solve then becomes 

W{a, XL)W{a, -Xr) = W{a, -XL)W{a, xr) . (23) 

By utilizing the asymptotic expansion of W{a,x) and W{a, —x) given in we can 
reduce eq. to the form 

^ sin sin = cos v^L cos , (24) 



up to terms of order h. In (j2 



1 o , . , TT 1^ 

(25) 



ipL = - aln(-a;L) + - + -$ 



and 



1 2 , 
'^R = -^^R - alnxi? + ^ + 2*^2 



k = Vl + e2™ - e™ (26) 
$2 = argr(^ + za). (27) 
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A trigonometric identity brings 



into the form 

2™ COs[ipR + ipi) 



cos(v3i? - v^l) = 
which upon using the definitions and collapses to 

cos C = vT+fr^™' cos 2lpo , 

where 

TT 



2(fo = ip^ + ip^ = -Xg - 2alnxo + - + $2 , 



(28) 



(29) 



(30) 



and Xq = ^/2/h, and we need only retain terms up to that order in h. Eq. 
determines the dependence of the eigenvalues a of (PT| on C. Solutions of this equation 
come in pairs, 



=F- arccos 



cosC 



, , (31) 

for each integer n. By explicitly constructing the wavefunctions associated with these 
eigenvalues, one can check that the solutions a^/a^ correspond to even/odd modes, 
respectively, where we remind the reader that we qualify modes as even or odd de- 
pending on their parity under the combined refiection {x,C) —>■ {—x,—C). We also 
recall that we can use ()31|) to compute the density of states of the inverted harmonic 
oscillator. Summing even and odd modes, the density of pairs is given up to terms of 
order 1/| ln^| in the double-scaling limit by 



p{a) 



dn 



da 



TC 



(32) 



In the limit, p diverges as | ln/i|. 

To proceed, we parameterize the solutions of as 

a± = ao + -[ai ± a2) , 



(33) 



where we define oq by the property fo{ao) = rni, i.e., is independent of C, and ai 
and 02 vanish in the double-scaling limit as 1/| In/ip and 1/| ln/i|, respectively, see eq. 
(IHUj). We can then expand for fixed ao to find 

1 , 



V9o(ai ± 02) 
where everything is a function of oq and 



TV - -V 02 



Via) 



arccos 



cosC 



.VTT 



-2na 



(34) 



(35) 
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Solving ()34j) yields the leading corrections to the eigenvalues in the double-scaling limit, 



-t;(ao) 

1 , 
--V a2 . 



(36) 



We have plotted this dependence of the eigenvalues (j33|l on C in Fig. El on page Cni 
The plot shows —a{C) for the first few eigenvalues closest to the top of the potential. 

We are now in a position to evaluate the finite contribution to the ground state 
energy of our matrix model that depends on the zero mode of the RR scalar C. As 
we have explained in the previous section, we turn on RR flux by filling even modes 
up to /i+ = /i + Q and odd modes up to /i_ = /i — Q from the top of the potential, 
similarly to Fig. |21 on page |H1 The ground state energy is then given by summing over 
all the eigenvalues (jHTj) in the parameterization (jHH)) (recall that in the conventions of 
this section, a is counted from the top of the potential downwards) 



-a2 



(37) 



Y,a+ + ^a^ = 5^(2ao + ai) + - (^0 + -ai) " 3' 

The rationale behind this particular splitting is the following. The first term in (j37j) 
is independent of the flux. It has the usual divergences computed in fSSJ. The second 
term depends on the flux, but is independent of C. It has a divergence which represents 
the bulk contribution of the background flux (jS)). We will analyze these contributions 
in more detail in the next section. The last term in ^7\i is the one of present interest. 
By plugging in the various definitions and using the density of states (j32p . we find that 
this term is finite in the double-scaling limit and given by 



^(C) = - / da^ p^ao) -03(00) = / da V9d(a)a2(a) 



1 

2^ 



— / da arccos 



cos C 



-2iva 



(38) 



which is the main result of our paper. We show a plot of this function in Fig. |3] 

The term (|38|) is the only C-dependence that survives the double-scaling limit. By 
curiosity, one may also ask for the next subleading term. We can easily extract this 
from the first term in (jHTj) and find it to be 



A'K\\nh\ V 



arccos 



cos C 



-2'jTa 



(39) 



15 



V{C) 




Figure 4: The ground state energy of the matrix model depends on the boundary condition 
asymmetry C. On the left, we plot the contribution (|38() which is finite in the double-scaling 
limit, for /i = .4 and Q = .1. On the right, we show the subleading term (|39|1 which goes to 
zero as 1/| \nh\. 

In the perturbative limit fi ^ oo, fjHH|l obviously reduces to the result (fT^ found 
in the previous section. The only exceptions are the regions of order e~'^^ around the 
minimum and maximum of V{C), where tunneling becomes dominant. It is curious 
to note that these non-perturbative effects actually lead to a large second derivative 
V" ~ e'"'^ at the minimum of the potential, in other words, the "mass" of the C- 
field appears to be of order e^^^^l We strongly emphasize, however, that the quantity 
V{C) we have computed cannot strictly be thought of as a mass term for C. In the 
spacetime interpretation, this potential energy is finite only after integrating over space, 
see eq. (fT^. Thus, the instanton contribution to the energy density vanishes, as could 
have been expected based on the fact that the instantons in question are D-branes 
localized at ^ oo in the Liouville picture ^HlEl- On the other hand, we note from 
the action Q that the kinetic term for the C-field goes to zero at — > oo, so that 
even an infinitesimal potential can conceivably lead to a dynamical stabilization of the 
constant part of C. More precisely, if we assume that the potential indeed comes from 
the instantons localized at infinite 0, we can imagine cutting out a finite part of the 
Liouville direction at large where the kinetic term for C is very small. In this region, 
there would be a finite energy density from V{C), which would stabilize the zero mode 
of C (its value at the cutoff) at the minimum of the potential. 

Another interesting feature of V{C) is that its minimum is at C = or C = vr 
depending on the sign of Q. This seems to be a property of non-perturbative potentials 
for axion-like fields also in higher dimensions, as for example the ones considered in [2]. 

The term (|39|) is also interesting. First of all, its minimum is at C = tt, which is 



16 



somewhat unexpected. (For symmetry reasons, the minimum can only be at C = 
or C = TT.) Moreover, ahhough it vanishes in the double-scahng hmit, one can check 
that its second derivative with respect to C is in fact infinite at the minimum C = ir. 
Thus, even though V{C) itself vanishes at Q = 0, this subleading term seems to fix C 
with high precision at C = vr also in the absence of fiux. This is true in particular if 
we imagine working with a cutoff Liouville direction. 

We have here computed the exact answer for the potential from the matrix model 
point of view. It would naturally be extremely interesting to understand some of the 
features of the potential from the spacetime perspective. 

Throughout the discussion, we have assumed that both Fermi levels are well below 
the top of the potential, i.e., /i > 0, and \Q\ <^ Only then is C a perturbative 
symmetry of the matrix model. When continuing /x to negative values (but keeping 
\Q\ <^ |yu|), Q should be interpreted as the magnetic flux ©, or alternatively as electric 
flux of the dual scalar C. Our formulas for V{C) still make sense in that case, but the 
perturbative symmetry is the one associated with C. We compute the potential V{C) 
in appendix 1X1 



4 Finite temperature partition function and dualities 

We will here follow the usual route developed in |^ for extracting flnite thermody- 
namic quantities from the double-scaled matrix model. In the matrix model (at flnite 
A^), we may consider the partition function 

Z = Tt e-^^^-'^^)) , (40) 

where f3 = 1/T is the inverse temperature and /i is the chemical potential. We note 
that we are here using rescaled variables that stay flxed in the double-scaling limit. 
Thus, H in (|4(Jj) is the Schrodinger operator (|2H). and /i > is measured from the top 
of the potential as in ()20|) . We will parameterize (3 = 271 R, as is appropriate when 
thinking of ()40|) as the path-integral with Euclidean time direction compactifled on a 
circle of radius R. 

An important quantity for deflning the thermodynamic limit of interest for string 
theory is the mean particle number 

A = (AT) = - Tr iVe-^(^-'^^) = ^ , (41) 
Z Oji 
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where F = \nZ/(3 is the grand-canonical free energy. More precisely, the double- 
scaling limit which defines non-critical string theory involves ^ — > 0, A ^ oo, keeping 
the chemical potential fi fixed. Since we are dealing with free fermions, it is easy to 
write A in the thermodynamic limit as 

/oo 
dapia)f{a), (42) 
-oo 

where p is the density of states and / is the usual Fermi distribution function 

/W=/K/^)= l + e^L(.-a) - (43) 



The diverging expression (j42|) for A can be studied more conveniently by taking 
two derivatives with respect to p. Using that df /dp = —df /da and f'{a) — * for 
a ±oo, one finds 

a^A f 

'da p'{a)f{a) . (44) 



We briefiy review the case of the bosonic string. Using (jH^ and the integral repre- 
sentation of the digamma function {r is often known as ip, but we have already used 
that letter), 

1 1 1 f°° t/2 

p'(a) = -Imr'(- + .a) = -y^ *-^smat, (45) 

as well as 

/'(a) = — ^ , (46) 

2 cosh 7!-R{p — a) 

we can write "Fourier integral" 

92 A 1 p , t/2Rp t/2p 

■ L / ■ i I sint. (47) 



dp^ 2ttp Jq smht/2Rp sinht/2/i 

For comparison with jOj, we note that we are here working with units in which a' = 1, as 
follows from the fact that the curvature of the matrix model potential is V = — 1, 
which is identified with the mass of the tachyon = —1/a'. The result fl47|) is 
invariant under T-duality, 

R^R' = ^, p^p' = ^p. (48) 

More precisely, recalling that A is related to the partition function via three derivatives 
with respect to p and a factor of i?, we see that 

Z{R,p) = Z{R!,p') (49) 
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holds for the bosonic matrix modeL 

When studying the OB theory, we may introduce two different chemical potentials 
/i+ and /i_ for even and odd modes, respectively. As before, /i± = fi ± Q, where 
Q is related to the RR flux background as we have explained in section |21 For the 
computation of the partition function, one has to be slightly careful to split the various 
contributions in the right way. Starting at finite N, we have 

A = ^(/(a+,/i+) + /(a_,/i_)), (50) 

where the sum is over all eigenvalues which are given by the solutions of (j^^ . One 
way to evaluate ()50p is to work with different "densities of states" for even and odd 
modes, defined by differentiating (jHlll . On the other hand, there is an even mode for 
every odd mode, so that the two densities should actually be the same. The answer is 
identical, but we find it cleaner to proceed by using the parameterization in ((221) • We 
obtain 

A = 5]^(/(ao,/i+) + /(ao,/i-)) + (/'(oq, /i+) -/'(ao,/^-)) 

f 1 (51) 
cia p(a)(/(a,/i+) + /(a,/i_)) + / da —v{a)[f'{a,^+)-f\a,^J)), 

where we use the same density of states p as in the bosonic case ()H2|1 . We can now 
again apply two derivatives to ()5H) . and use the representation. 



1 , TT sinhvra 1 , t 2 

—v"ia) = 5 — = — / dt — ^— sm at . (52) 

27r ^ ^ 4cosh2 7ra 2Tr Jq cosht/2 ^ ^ 



In this expression, we have set C to the minimum of its potential (O/vr depending on 
the sign of Q). It is possible to compute the Fourier transform of v'{a) also for general 
C, but we relegate the somewhat cumbersome expressions to the appendix. We now 
find 



d^A 1 p t/2R 



-Im (e*^+* + e'"-') + , Im (e*^+* - e*'^"*) 



9/^2 27r Jq sinht/2R [sinht/2 ^ ^ cosht/2 



1 , t/2R 
' dt ' 



Tc Jq sinht/2i? 



t/2 . ^ t/2 . ^ ■ 

sm fit cos Qt H — — cos fit sm Qt 



sinh t/2 cosh t/2 



(53) 



For Q = 0, and remembering that a' = 1/2 in our conventions for the supersymmetric 
model, this reduces to the expression given in j^. The reader might wonder about the 
second term in the square brackets in ((^. which appears to be odd under Q —Q, 
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whereas the partition function should naively not depend on the sign of the flux. The 
resolution of this puzzle is that we have not included the dependence on the zero mode 
of C, the minimum of whose potential switches from C = to C = vr as we change the 
sign of Q. Taking this into account makes symmetric under Q ^ —Q. (See eq. 
in the appendix.) 

The result fj53|) is the complete non-perturbative result for the partition function 
as a function of /i, Q (and C). One can think of the first term in the square brackets 
as the non-perturbative resummation of the asymptotic series in that is defined 
in string perturbation theory. In fact, this asymptotic series is essentially the sum of 
the two series 

a;^,,,(/x, q) = a;^Ai^^) + AiAf^_) , (54) 

where the perturbative expansion of Apgj,^(/i) is identical to the one of the bosonic 
string 



A''(n) = — I dt ^/^^ „ Ime''^* 



1 



smht/2R sinht/2 

_i ^ 1 I R I I R 

Ji y J (2m)! y |22'=-2| |2^-^^-2| '7,^f '^^'"-^^ 
27r ^^/x2™+i(4i?)'™^ {2k)\{2m-2ky. 

(55) 

where the are the Bernoulli numbers. On the other hand, one can see that the 
second term in square brackets in (j53|l is analytic in fi and its asymptotic expansion in 
vanishes, thus indicating its non-perturbative nature. To make this explicit, we 

write 

A:^if^,Q) = Alifi+)-A:^ifi.) (56) 

with 

The essential feature of this integral is that the integrand is even under t —>■ —t. We 
can then close the integration contour in the upper/lower half plane (depending on the 
sign of /i) and write the integral as a sum over poles. Explicitly we have for yU > 0, 

A" (,,\ _ V /r_lV"+l ^^^^^^ -27rnMi? , Nn+l - W^^'^ / R 2^(n-i)/.l 

cosnvr/? ^ sin{n - \)ti / R J' 

(58) 

and a similar expression for yU < 0. We see from (j58|) that there are two types of 
non-perturbative contributions. One type comes from the poles of the ^^^^^^^ factor in 
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()57|) and begins at order e~'^^.^ The spacetime interpretation of these contributions 
are D-instantons, i.e., ZZ-brane jlS] in the Liouville direction together with Dirichlet 
boundary condition in the Euchdean time direction. The other type of contributions 
comes from the poles of the thermal factor in (|H7jl and begins at order e"^''^''. Ob- 
viously, the spacetime interpretation of this term is a DO-brane winding around the 
compact Euclidean time direction, i.e., ZZ-brane together with Neumann boundary 
condition on the circle. It is interesting to note that something special happens to 
the structure of the non-perturbative series when the radius R is rational. In that 
situation, the two types of poles collide and combine into a single contribution (this 
limit is built into ()58|) automatically). It would be nice to understand this structure 
from the Liouville point of view. The relation between ZZ-branes and non-perturbative 
contributions in the bosonic matrix model has also been discussed in ^\ . 

We stress that it is really sensible to write an instanton expansion for this second 
term in because it is odd under Q — > —Q and hence vanishes at Q = 0. The 
corresponding amplitudes vanish in perturbation theory, and receive contributions only 
from the instantons. This property is in distinction to the first term in (jSHI), whose 
perturbative expansion is non-zero. From this asymptotic expansion, one can merely 
reconstruct that its non-perturbative ambiguities set in at order g"^'^'^^ and e"^'^^, but 
it is not sensible to write an instanton expansion for these ambiguities. 

With these results in hand, we now proceed to study various dualities that we 
expect our theory to enjoy. We begin with T-duality. The partition function of type 
OA theory has been computed in [H] using the matrix model. It is found to be 

a^A 1 P , t t/2R . I 

= — dt — ; ^- — sinute-l^l*, 59 

d/j,^ 2tt Jo sinht smht/2R ^ ' ^ ^ 

where q is the integer flux parameter of the OA theory measuring the net number of 
DO-branes. Is is easy to see that the results for q = Q = are related by the simple 
T-duality 

Za{Ra, fJ'A, q = 0) = Zb{Rb, /is, Q = 0) , (60) 
where the parameters are related by 

Ra = 77^ = ^ ' f^A = 2RBfiB = -^== f^B ■ (61) 

2Kb Kb \/ol 12 

It is worthwhile to point out that (j61|) is a slightly non-standard transformation rule 
on the string coupling, which in the usual type II context transforms as in (pSj) . It 

^Interestingly, for C ^ 0, these contributions start at order e"^'^'^, see appendix IbI 
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would be interesting to understand why the string couphng has to be normahzed in a 
different way in this two-dimensional type context. 

If we want to check T-duality also for non-zero flux, the form of the above expres- 
sions suggests that we should try to Wick rotate q iQ and supplement (jF)T|) with 

iq = QA = QB^ = ^- (62) 

The fact that the flux transforms with a factor of i is quite suggestive of a timelike 
T-duality. Namely, recall that the flux-parameter Q in type OB is for an electric flux, 
which is a timelike one-form. Rotating to Euclidean time, the flux becomes imaginary. 
T-duality turns the one-form into a "zero-form flux", which is really just the dual of 
the 2-form flux of OA. But now undoing the analytical continuation does not give back 
the factor of z, so that the flux stays imaginary. We will then take (jUT|) and ()62|) as 
a hypothesis for T-duality, and compute the partition function for the OA model at 
imaginary q = iQ. We will be able to match the perturbative expansions of the two 
theories, but we will not be able to find agreement at the non-perturbative level. 

In the type OA matrix model introduced in jH], the Schrodinger equation for the 
eigenvalues (the analogue of (EH) is (see also jlSllllI for early studies of this matrix 
model potential) 

/ 1 x'^ q^ \ 

\-dxxdx + - r - a)^(a,x) = . (63) 

Writing z = ^ and iIj = z^/^e~^/^/, we obtain the confluent hypergeometric equation 

(z9,2 + (l + g-.)9.-(i + |-|))/ = 0, (64) 

with well-known solutions. The problem with the analytical continuation q ^ iQ 
comes from the boundary condition at x = 0. Indeed, the two solutions of ()63j) behave 
as ip ^ x°' with = near x = 0. We see that for > 1, only one of the solutions 
is normalizable, Jg x^"^^ < oo and there is no need for a boundary condition. For 
< 1, the singularity is of limit-circle type and we need an extra boundary condition 
to make sure that the differential operator (jUHj) is self-adjoint. The usual procedure 
for = is to require that the Laplace operator in the plane be self-adjoint. But for 
g^ = —Q^ < 0, the solutions are actually oscillatory near x = 0. The possible real 
boundary conditions come in a one-parameter family, 

^ ^ e'^/^x'^ + e-'^^^x-''^ as X ~ 0, with 6 real. (65) 
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Using the traditional Kummer notation for the confluent hypergeometric function, we 
can write the solution of ()63|) satisfying the boundary conditions ()65|) as 

ilj{a,x) = e'^/^x'^e-'/^<l>{^ + f-f,l+iQ]z) + e-'^/'^x-'^e-'/^<!>{^-f-f,l-tQ;z) . 

(66) 

By working through the asymptotics of the $ function, one finds that this solution 
behaves near x —>■ oo as 

(67) 



1 OC TT 

'ipia, x) ~ const. — cos( a In —= ip) 

X M V 2 4 ^ 



where the phase shift (f is given by 



Lf = ar^ 



gi(5+Qln2)/2^X^^^^^^)_g~7rQ/4 _^ ^-»(3+Q In 2)/2 T{l-lQ) ^^q/^ 



r(i + f + f) 



r(i-f + f) 



After absorbing an a-independent phase into 6, this reduces to 



if = arg 



gj5/2g-7rQ/4 g-i5/2g7rQ/4 
+ 



r(i + f + f) r(i-f + f) 



(68) 



(69) 



It seems hard to find an explicit integral representation for the density of states 
similar to (jl^ . But noticing that 



r(i-f + f) 



r(i + f + f) 



e™ + e 



-ttQ 



(70) 



we can expand the problem for large /i ^ 1. We write (f = (fi + (f2 the sum of two 
terms 



-2 [argr(| + f + f)+argr(l-f + f)] 



{p2 = arg[w + 



(71) 



where w = [e'^ e'^'^/^T f + f) /T{l + f + f)] ^''^ has absolute value exponentially 
close to one. The first term in (f7T|) gives after the usual steps 



1 , , 



t/2R 



■ sin /it cos - 



(72) 



sinht sinht/2i? 

which is precisely the T-dual of the first term in (j^Hj) . as we have claimed. This 
fact that the asymptotic expansions of the two partition functions agree to all order in 
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perturbation theory can also be seen by analytically continuing directly the result ()59|) . 
Indeed, as we have seen above, only the even part of e"''* will lead to a perturbative 
series. Thus, the perturbative expansion can be obtained by replacing the exponential 
with coshgt, which leads to a convergent integral for sufficiently small q. Now rotating 
q ^ iQ directly shows the agreement with the first term in ()53|) . But of course, it is 
safer to first analytically continue the problem as we have done above. 

What about the non-perturbative contributions? On the OA side, these come from 
the second term in We find for large a ^ Q 

(/P2 ~ e-^"sinh7rg tan(f - f Ina) , (73) 

where we have again absorbed all a- independent phases into 6. We see that the first 
non-perturbative contribution to the OA partition function at imaginary fiux Q = iq 
is of order ~ e~'^'^ sinh ttQ and depends in a somewhat strange way on the boundary 
condition 6. This result can indeed be matched with the order of the first correction 
in (j58p after using the T-duality transformation rules (jblj) . The fact that we need the 
additional boundary condition at x = is suggestive of trying to match it with C. 
We have computed the corresponding dependence of the first instanton contribution 
on C from the formulas in the appendix but have not been able to match this with the 
dependence of (ffH|l on 6. And the prefactors of the exponentials do not show an exact 
match either. 

Moreover, one has to ask for the presence of the second type of non-perturbative 
contributions, which begin at order e~^^ on the OB side. On the OA side, this would 
be a contribution of order e~'^^^ and has to come from the thermal density factor 
1/(1 + e27rfi(M-a)~)^ these contributions start at order e"'^'^^^, in disagreement with 
expectations. It therefore seems unlikely to us that the naive prescriptions that we 
have followed above can reveal an exact T-duality between OA and OB, once the fiux 
is turned on and non-perturbative effects are accounted for. 

Non-perturbative violations of T-duality have been noticed before, as for example 
in jlTj. In this paper, the duality group of the heterotic string on K3xT^ is studied. 
It is argued that the failure of naive T-duality at the non-perturbative level should be 
viewed not as a breakdown of perturbative intuition, but rather as a deformation of 
the latter. In the present context, this would suggest to deform the T-duality relations 
f)61|) by terms of order e~^^ and e~'^^^ such that the partition functions agree exactly. 
It would be worthwhile to study this proposal further. 
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Before leaving the subject, we also wish to clarify what we mean by "the type 
OA matrix model at imaginary g" . The definition of the OA matrix model in ^ uses 
rectangular (A^ + g) x A^ dimensional complex matrices, and a U{N + q) x U{N) 
gauge group which removes the phase of the "eigenvalues" . The holographic intuition 
behind this definition is that the matrix model is the worldvolume theory of a large 
number of -DO and DO branes, and q is the net Z)0-brane charge. Leaving aside the 
fact that this theory is rather unphysical to study at imaginary q, it is not even clear 
how to define the matrix model mathematically! What we have in mind here is an 
alternative definition of the OA matrix model, in which the q leftover DO-branes have 
been dissolved into background flux. Similarly to the discussion at the end of section 
12 the tachyon T on a general DpDp system, which is a complex field, couples to the 
RR C^P~^^ form potential such that a vortex of T carries one unit of lower-dimensional 
charge. This corresponds to a term § daxgT A C*^^"^-* = / d{da,rgT) A C^^~^^ in the 
worldvolume action. Integrating by parts, this becomes /rfargT A F^p\ Thus, in the 
type OA theory in two dimensions, if we turn on two-form flux F^'^^ = *F^^^ oc q, we 
obtain a term q J J^(argA) in the quantum mechanics of the complex eigenvalues A. It 
is easy to see that after gauge fixing, such a term corresponds precisely to a shift of 
the angular momentum in the eigenvalue plane by q units, just as in (jUHj) . This way of 
defining the OA model does not pose any obvious obstacles to g ^ iQ, except for the 
fact that it might be difficult to define the model at finite A^ because the Hamiltonian 
is unbounded from below. 

Another interesting duality jH] of the type OB theory in two dimensions is "S- 
duality" fi —>■ —fi and C — > C, as we have mentioned several times already. Our results 
for the partition function at C = are obviously invariant under /i —fi, Q —Q 
which is simply particle-hole duality of the matrix model eigenvalues. It is somewhat 
less clear what becomes of our non-perturbative results involving the zero mode of 
C. While the formulas do not suffer any obvious pathology, C as described in section 
12] is not a perturbative symmetry for fermions above the top of the potential. For 
such states, the perturbative energy levels are left/right moving, and non-perturbative 
effects are reflections off the potential. The perturbative symmetry associated with 
these energy levels is precisely the zero mode of the dual variable C, and one can check 
that the non-perturbative dependence of eigenvalues on C is precisely as in (jHU, with 
a —a. This is discussed in detail in appendix 1X1 

We can also notice yet another duality from our matrix model computations. The 
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partition function (which is the third integral of (jSHj) with respect to /i) is also invariant 
under yU„ — Z^-, ^-c, applying particle-hole duality to the odd modes only. In terms 
of and Q, this corresponds to exchanging NS and RR backgrounds fi ^ Q, a rather 
intriguing duality indeed! The T-dual of this type of duality has been noticed before 
in the OA context in |8j. It is much more manifest in the OB model. 

To summarize, the OB matrix model depends on two parameters fi and Q, or equiv- 
alently /x+ and but models which differ only by a sign of fi± are equivalent to each 
other. If \Q\ <^ the matrix model is the non-perturbative definition of a OB non- 
critical string theory. Depending on the sign of /i, Q is interpreted as either electric or 
magnetic flux with corresponding perturbative gauge symmetry C or C. It would be 
interesting to understand whether there is also a perturbative string theory description 
in the case <C \Q\. An analog of the perturbative gauge symmetry does not seem 
to exist in that case. 

5 Bosonization and amplitudes in RR flux background 

According to our identification in section the matrix model with different Fermi 
levels for even and odd modes provides a dual description of type OB string theory 
in the RR flux background (jSj). This formulation is extremely simple and makes it 
possible to compute scattering amplitudes and a complete S-matrix for this string 
theory. In this section, we will describe the formalism for studying these amplitudes, 
leaving a detailed study of their properties for future work. Our approach follows the 
older works [101 1^ on the bosonic string and the more recent paper [7j on the type 
OB model (with vanishing RR flux). The use of matrix models to study non-critical 
strings in flux backgrounds was first proposed in the context of type OA in jH]. 

In relating matrix quantum mechanics to spacetime physics in two dimensions, one 
identifies collective excitations of the Fermi sea of matrix eigenvalues with the bosonic 
fields of the corresponding string theory. This identification includes a non-local field 
redefinition which is implemented by the so-called "leg-pole factors" in the scattering 
amplitudes, but we can ignore this subtlety for our purposes. For the OB model, it 
was proposed in j2J|H] to identify the even fluctuations of the Fermi surface with the 
tachyon field T and odd fluctuations with the RR scalar C. In j7j, it was pointed 
out that unitarity of the S-matrix requires the inclusion of "solitonic sectors" into the 
Hilbert space of asymptotic states. These solitonic sectors are created by two additional 
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bosonic fields, which we will call S'^ and S~ , which carry C-field charge, and have no 
perturbative spacetime interpretation. 

More precisely, as in PU], one begins by introducing second quantized fermion fields 



oo 

iwt 



J — oo 

where the ip'^{uj, A) are a complete basis of solutions of the classical field equation ()18p or 
fl21|) . and where we identify uj = —a to conform with other conventions in the literature. 
For fixed u, different such bases are appropriate in different situations. Of primary in- 
terest for defining scattering amplitudes are in/out bases which are defined by requiring 
specific asymptotic behavior. We will denote solutions which are in/outgoing on the 
left and right of the potential by subscripts L and R, respectively. The 'in' and 'out' 
bases are related, as usual, by the S-matrix, 



^out 




(75) 



where we have chosen the phases of L and R modes such that the Z2 symmetry of the 
potential is manifest {i.e., S commutes with (iq))- textbook methods, one finds 
the refiection and transmission coefficients to be 

R{iu) = -Le-^'^/2r(^l _ ^ y^^) _ ^e"^i?(cJ) . (76) 
V 27r 2 



We stress that ()74p is the full "interacting field", where in the present case the only 
interactions are with the external potential, and does not depend on the conventional 
choice of basis. For example, for the purpose of defining the vacuum, it will be manda- 
tory to introduce other bases of modes with fixed parity. 

But let us for the moment focus on the asymptotic left /right bases. Both in and 
out modes yield the usual algebra of creation and annihilation operators, 

{b^iu),biiLu')} = 5iuj-uj') 

{bn{u),bl{uj')} = 6{u;-u;') (77) 
and {bf (cu) , bf {to')} = for all other combinations. 

To complete the definition of the theory, we have to fix a representation of this algebra, 
as we will do a little later. But independent of this representation, we can introduce 
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bosonic operators as fermion bilinears such as 

/oo 
-oo 

/oo 
biiObni^ + uj) 
-oo 

As usual, these bosonic operators satisfy the relations of a current algebra^ 

(79) 

[aRR(t^i), aRR(u;2)] = ujiS{uji + ^2) ■ 

In jHlini, it was proposed to identify the spacetime tachyon field T with the even 
combination oll + Qrr and the RR scalar with the odd combination orr — oll- But 
as pointed out in jT,, the bosonic S-matrix computed with these fields cannot possibly 
be unitary. Indeed, it is easy to see that we can define two further fermion bilinears. 



/oo 
-00 

/oo 
bi{ObL{i + 
-00 



00) 



^0) 



which complete (|79|) to a f/ (2) current algebra 

[aRR(t^i), aRL(u;2)] = aRhi^Ji + UJ2) 
[aLL(^i),aRL(u;2)] = -aRL(t^i +^^2) 

(81) 

[aRL(t^l),aRL(t^2)] = 

[aRL(cJi), aLR(u;2)] = aRR(t^i + 002) - au^iyji + LO2) + uji5(ui + U2) , etc. 

In this algebra, the tachyon T ~ oll + orr is one of the f/(l)'s and the RR scalar 
C ~ flRR — oll together with the "soliton operators" ~ olr and S~ ~ orl generate 
an SU (2) algebra. The main point of [7j is that because the interactions (scattering 
of fermions off the potential) do not respect this symmetry algebra, it is inconsistent 
to restrict the Hilbert space of asymptotic states to a fixed charge under the U (1)^ 
currents (ffn|) . 

Going back to the fermionic picture, we recall that we still have to specify the 

vacuum of our theory. This amounts to splitting the fermionic algebra (f77j) into creation 

^Actually, obtaining this algebra requires appropriately regulating in the UV, as is the case in our 
situation. 
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and annihilation operators, in other words, the choice of a Fermi leveL As we have 
explained in detail in section |2l the only non-perturbatively stable possibility is to 
choose separately Fermi levels for the even and odd modes, related to the L/R modes 
by 

&±H = ^(&lM±&rM). (82) 

Thus, we choose /i+ = /i + Q, /x- = /i — Q, and assume for simplicity that /i+ < and 
11- < 0. The vacuum |0) is defined by^ 

7^ ' (83) 
b\{uj)\Q) = Q uj<ij±. 

As we have mentioned, for fi^ = fi_, the bosons defined in ()78|) have the interpretation 
of collective excitations of the Fermi surfaces which are coming in from the left or from 
the right, while the off-diagonal bosons (fHOjl create solitonic sectors. For Q 7^ 0, we 
loose a good semiclassical picture of the ground state, so this interpretation is not so 
obvious anymore. 

Finally, we write out explicitly the two-point functions computed using the defini- 
tions (f7S|). (|Hm) with respect to the vacuum defined in (jS^ . The formulae are simplest 
in terms of a±± = bl.b±, which can be related to T, C, S*^ and S*^ by working through 
the definitions (IH^. We find 



(a-|(a;i)aV\(cu2)) = 6{uj^ + u^) / d^i R+{^i + uJi)Rl{^i) 
(a°"l(a;i)ai_\(^2)) = 5{cu, + ^2) P d^i R-{^i + ^i)i?;(ei) ^^^^ 



(a7!(a7i)a-_(c^2)) = 0, etc., 
where the are the eigenvalues of the free fermion S-matrix ()75|) 

R± = R±T = - luj) (^ze-""/2 ^ j (§5) 

6 Open Problems 

In this paper, we have studied aspects of the RR sector of type OB string theory in two 
dimensions, using the matrix model as a non-perturbative definition. The crucial step 



^In the conventions of this section, /i < when the Fermi levels are below the top of the potential. 
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was the precise identification, following [B|, of the RR flux background as the difference 
of Fermi levels for even and odd modes of the inverted harmonic oscillator potential 
of the double-scaled matrix model. We have also identified the zero mode of the RR 
scalar C in the finite- iV matrix model as the mismatch of the perturbative (in K) energy 
levels on the two sides of the potential. 

We have then given three applications of this identification of RR fields in the 
matrix model. Firstly, we have computed the potential energy for the zero mode of 
C, which is induced by turning on flux, and stabilized by non-perturbative effects. 
Our result is the complete answer to this question, the analogue of which in higher 
dimensions is of great interest for the construction of realistic string theory vacua. 
Secondly, we have computed the finite-temperature partition function of the OB string 
in the flux background. We have seen that T-duality with the OA model is realized at 
the perturbative level, but appears to be violated by non-perturbative effects. Thirdly, 
we have outlined how one can use the matrix model to study scattering amplitudes 
of OB strings in RR backgrounds, following |inj|lTl|Hl[7j . Understanding strings in RR 
backgrounds is also an important problem in higher-dimensional string theories. 

To conclude this paper, we list here a few open questions which have been raised 
by our results. 

■ Quite obviously, it would be extremely interesting to understand the spacetime or 
string theory computation of our result for the potential V{C). There, the origin of 
the potential are non-perturbative effects associated with D-instantons, and it would 
be interesting to see how they contribute to V{C). 

■ More generally, it would be interesting to compute some other terms in the effective 
action and to see the structure of non-perturbative effects, which could in particular 
provide guidance for similar questions in higher dimensions. 

■ It would be nice to understand better what happens to T-duality with OA for non- 
vanishing flux. 

■ The matrix model gives complete control over interpolations between different dual 
perturbative string theories. It would be worthwhile to understand what happens at 
the transition region, e.g., when /i_ goes through zero at finite string coupling. 

■ What is the analytical behavior of scattering amplitudes in non-zero RR background? 
Can one derive the identification of RR fields directly from a discretization of super- 
Riemann surfaces? 

■ According to ^H], the old-fashioned c = 1 string at the self-dual radius is related to 
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the topological string on the conifold. It would be interesting to see whether our model 
with unequal chemical potentials corresponding to non-zero RR background could have 
a similar interpretation in the context of topological strings. 

Needless to say, we hope to come back to some of these problems in the future. 
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Appendix 



A V{C) = V{C) 

This discussion is mutatis mutandis almost identical to the one given in section |21 and 
El When the chemical potential fi is above the top of the matrix model potential, the 
perturbative energy levels are left /right-moving moving waves instead of localized on 
one side of the potential. The magnetic gauge symmetry associated with C shifts the 
phase of left movers by e**" and the phase of right movers by e"**" . We can implement 
C in the finite-A matrix model by imposing periodic boundary conditions on the 
wavefunctions with a phase shift of e**". This means that we study the Schrodinger 
equation and impose 

6*^/^(0, xo) = e~*^/^^(a, -xo) 



86) 



where xq 
equation 



— i> 00 in the double-scaling limit. An alternative is to study the 

C ^2 



Ox — « H 

2xoJ 4 



ip{a, x) = . 



7) 



with strictly periodic boundary conditions on ip. This description is of course gauge 
equivalent to ()8(i|l . but makes it more manifest that left /right movers pick up a phase 
shift when traversing the potential from — xq to Xq. This implementation of C is 
appropriate if we construct the matrix model using unitary matrices with periodic 
potential. 
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Upon writing the general solution of ()21|) in terms of parabolic cylinder functions 
W{a,x) and W{a, —x), eq. (jH^ becomes the eigenvalue equation 



e'^/^W{a, x) - e-'^/^W{a, -x) e'^/^W{a, -x) - e-'^/^W{a, x) 

e'^/^d^W{a, x) + e-'^/'^d^ {W{a, -x)) e'^/^d^ {W{a, -x)) + e''^/^d^W{a, x) 

Using the asymptotics for x ^ \a\, 

[2k 

W{a,x) ^ J — cos if W{a, —x) ^ y — simp (89) 
with (f = — alnx + ^ + |$2 and k = ^/l + e^™ — e^", we obtain from 



cos C = vT+e^™ sin 2ipo , (90) 

where (y^o = for x = xq = \/2/h. Up to a shift of (po by it /A, this is exactly the 
"S-dual" of and shows that the non-perturbative potential one would compute 
for the magnetic variable C is identical to the one we found for the electric variable 
C. Which potential is sensible to compute depends on the sign of fi, i.e., whether we 
interpret Q a electric or magnetic flux from the point of view of perturbative string 
theory. 



B y{C) at finite temperature 

In this appendix, we compute the Fourier transform of v'{a) in in order to obtain 
a closed form expression for the non-perturbative piece of the OB partition function at 
finite temperature and for general C. The quantity of interest is 

/•oo ^ 

K^{fi) = J da—v'{a)f'{a,fi), (91) 

where 

2 cosh 7ri?(/x - a) 2itJ_^ smht/2R 

and 

1 n . cosC 

— V (a) = , 93 

2n ^ ' 1 + Vsin^C+e-2™ ^ ^ 



The integrand in 



00 



da e-*^^!! ^"^^ f94) 

1 + Vsin^ C + 
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has poles at e~^'^° = —1 and branch points at e~'^'^'^ = — sin^ C. It is convenient to 
place the cuts for n > from a = {tt- + — ^ Insin^ C to a = (n + + oo. We 
can then evaluate the integral for t > by closing the contour in the upper half 
plane, picking up the poles on the imaginary axis and integrating back and forth along 
the cuts. Paying attention to signs, we find the contribution from the poles to be 

-1 cosC cosC 



_ -(n+i)t 

„=o 2vr^ " Vsin^C-l |cosC|_; 



— -L r^„_LlU (JUS O UUS V^/ ,N 



(95) 
cosC 1/2 



I cosC| cosht/2 

To evaluate the contribution from the cuts, we substitute a = + |) ^ — In sin^ C — 
2^ Inz, with z G [0, 1], and are reduced to 

f'^ dz -2 iu/ . 2^\-T- -;zsin^C cosC 

/ (-1) e ^"+2-'*(sm^C) 2^ -_- , ^ 

n=QJo ^ l-zsin2C|sinC|yr^ 

= - |sinC|cosC(sin^C)~^ C dz z'^ il - zY^l^il - zsin^ C)-\ (96) 

TT cosht/2 ^ Jo 

We recognize the remaining integral as the particular hypergeometric function 

r(i - f )r(i) ^, ^t 3 It 

V 27ry V2y.^/^^ ;sin2C) 



r(l-f) ^' 2vr'2 2vr' 

= (2 " ^) Icl'sCllsinCI ^(^ ~ 1^' ^) ^-^^(2 ~ ^ ' ^) ^^^^ 
which defines the incomplete Beta function according to the usual sources. Pulling 
things together, we can finally write 



w / X 1 ^ , V2i? 1/2 cosC , , 

A' /i = —Re / dt e-'^^ — / , ^ — — 1 + g(t, C)] , 9^ 

"P^'^^ 27r io sinht/2/? cosht/2 cos , ;j , v 



where 



«(''-)4(^l)K-l'^)--^(i-|4)- (-) 

The structure of the non-perturbative instanton contributions that follows from these 
calculations is quite interesting. As we have seen in the main text the first 

instanton term for C = is of order e'^^. On the other hand, when C is non-zero, it 
can be seen from that non-perturbative effects set in at order e^^'^^ for fi —>■ 00. 
(But diverge as 1/ sinC for C 0.) This behavior is easy to understand in the matrix 
model. When C = 0, the perturbative energy levels on the left and right match, and 
the simple tunneling amplitude is of order e~'^^. When C 7^ 0, the perturbative levels 
do not match, and one has to tunnel forth and back to find a non-zero contribution. 
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C c=l at R = l 



In the bosonic c = 1 matrix model, many formulas simplify at the self-dual radius 
R = 1. The supersymmetric model is not mapped onto itself under R ^ as we 
have seen in the main text. But we can again evaluate some formulas more explicitly 
when we set -R = 1. For example, the perturbative part of the finite-temperature 
partition function is essentially the bosonic string result. 

2 



KM-lfdt^-^) sm,t 



— Re9;(/x^(l-^/i)) 



1 

leading to the well-known asymptotic expansion 



(100) 



2.a;;„.M~^ + I: '' (101) 



n=l 



In the present case, we of course have to add together two such terms with fi = fi±. 

We can also compute the non-perturbative part of the partition function exactly at 
R = 1. Consider (again for one of fi±) 



/•oo ^ 

^np(yu) = / da —v'{a)f{a,fi) 



oo 

°° , 1 cosC 

da 



which setting y = e'^'^^ and z = e~^'^" can be reduced to an elementary integral 



°° , 1 1 cosC 
dz 



27r Jo l + yz 1 + z Vsin^ C + z 
1 cosC 1 



2tc I cosC| y — I 



2 arcsin ( I sin C I y/y) — TT 
2 arcsin I sin C| — TT — ^/y I cos C| - 



Vl-l/sin^C 

(102) 

where an appropriate branch of the arcsin is understood once ysm^ C > 1. 
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